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The reduction operators, i.e., the operators of nonclassical (conditional) symmetry, of (1 + 1)- 
dimensional second order linear parabolic partial differential equations and all the possible reduc- 
tions of these equations to ordinary differential ones are exhaustively described. This problem 
proves to be equivalent, in some sense, to solving the initial equations. The "no-go" result is 
extended to the investigation of point transformations (admissible transformations, equivalence 
transformations, Lie symmetries) and Lie reductions of the determining equations for the nonclas- 
\& \ sical symmetries. Transformations linearizing the determining equations are obtained in the general 

case and under different additional constraints. A nontrivial example illustrating applications of 
reduction operators to finding exact solutions of equations from the class under consideration is 
p I ■ presented. An observed connection between reduction operators and Darboux transformations is 

discussed. 

j * , 

c3 . 1 Introduction 

s 

The notion of nonclassical symmetry (called also Q-conditional or, simply, conditional symmetry) 
[ was introduced in [2] by the example of the (1 + l)-dimensional linear heat equation and a 

; particular class of operators. A precise and rigorous definition was suggested later (see, e.g., [HI 

IHl HI]). In contrast to classical Lie symmetry, the system of determining equations on the 
j — , coefficients of conditional symmetry operators of the heat equation was found to be nonlinear 

and less overdetermined [2]. First this system was investigated in [42] in detail, where it was 
partially linearized and its Lie symmetries were found. The problem on conditional symmetries of 
the heat equation was completely solved in [7], see also [6]. Namely, the determining equations 
were obtained in both the cases arising under consideration and then studied from the Lie 
symmetry point of view and reduced to the initial equation with nonlocal transformations. 
The maximal Lie invariance algebras of the both sets of the determining equations appeared 
isomorphic to the maximal Lie invariance algebra of the initial equation. (Well later few of these 
results were re-obtained in [IS].) The results of [7J were extended in [22J, [23] to a class of 
linear transfer equations which generalize the heat equation. Thus, for these equations the "no- 
go" theorems on linearization of determining equations for coefficients of conditional symmetry 
operators to the initial equations were proved in detail and wide multi-parametric families of 
exact solutions were constructed with non-Lie reductions. It was observed in [43] that the 
proof of the theorem from [7] on reducibility of determining equations to initial ones in case 
of conditional symmetry operators with vanishing coefficients of dt are extended to the class of 
(1 + l)-dimensional evolution equations. This theorem was also generalized to multi-dimensional 
evolution equations [24J and even systems of such equations [30] . 

The conditional invariance of a differential equation with respect to an involutive family of 
/ vector fields is equivalent to that any ansatz associated with this family reduces the equation to 
a differential equation with the lesser by / number of independent variables [44] . That is why, we 
use the shorter and more natural term "reduction operators" instead of "operators of conditional 
symmetry" or "operators of nonclassical symmetry" and say that a family of operators reduces 
a differential equation in case of the equation is reduced by the associated ansatz. 
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In this paper we investigate the reduction operators of second-order linear parabolic partial 
differential equations in two independent variables, which have the general form 

Lu = u t — A(t, x)u xx — B(t, x)u x — C(t, x)u = 0, (1) 

where the coefficients A, B and C are (real) analytic functions of t and x, A ^ 0. These co- 
efficients form the entire tuple of arbitrary elements of class ([I]). We justify the partition of 
the sets of reduction operators into two subsets depending on vanishing or nonvanishing of the 
coefficients of cV Usually this point is missed in the literature on conditional symmetries. After 
factorization by the equivalence relation between reduction operators, we find the determining 
equations for the coefficients of operators from both the subsets. All the possible reductions 
of equations from class (pQ) to ordinary differential equations are described. Different kinds of 
"no-go" statements on reduction of study (including solution) of the determining equations to 
the corresponding initial ones are obtained for equations from class ([I]). In particular, the point 
transformations of all kinds in both the classes of determining equations (admissible transforma- 
tions, transformations from the associated equivalence groups, Lie symmetry transformations) 
are induced by the corresponding point transformations in class ([1]). Lie solutions of the deter- 
mining equations prove to admit nontrivial interpretations in terms of Lie invariance properties 
of the initial equations. An example on application of reduction operators is presented. It shows 
that in spite of the "no-go" statements nonclassical symmetry is an effective tool for finding 
exact solutions of partial differential equations. 

There are a number of motivations inducing us to carry out the above investigations. Class (pQ) 
contains important subclasses that are widely applied in different science (probability theory, 
physics, financial mathematics, biology, etc.). The most famous examples are the Kolmogorov 
equations (C = 0) and adjoint to them the Fokker-Planck equations (A xx — B x + C = 0) which 
form a basis for analytical methods in the investigation of continuous-time continuous-state 
Markov processes. (The other names are Kolmogorov backward and Kolmogorov forward equa- 
tions, respectively.) The first use of the Fokker-Planck equation was the statistical description 
of Brownian motion of a particle in a fluid. Fokker-Planck equations with different coefficients 
also describe the evolution of one-particle distribution functions of a dilute gas with long-range 
collisions, problems of diffusion in colloids, population genetics, stock markets, quantum chaos, 
etc. Due to their importance and relative simplicity, equations from class (pQ) are conventional 
objects for studies in the framework of group analysis of differential equations. Lie symmetries 
of these equations were classified by S. Lie [13] • The (I+I)-dimensional linear heat equation is 
often used as an illustrative example in textbooks on the subject |17| and a benchmark example 
for computer programs calculating symmetries of differential equations [TT]. It is the equation 
that is connected with the invention of nonclassical symmetries [2]. First discussions on weak 
symmetries also involved the linear heat equation and a Fokker-Planck equation [191 [37] . At the 
same time, all previous studies of nonclassical symmetries of equations (pQ) were not systematic. 
Only a few equations and single properties were considered. 

Results of [5j [7J [221 [27J are extended in the present paper mainly in two directions. Firstly, 
the entire class (pQ) is regularly investigated with the nonclassical symmetry point of view and, 
secondly, non-evident properties of point transformations and Lie reductions of the determining 
equations are found via involving admissible transformations in the framework of nonclassical 
symmetries. 

Our paper is organized as follows: Necessary notions and statements on nonclassical sym- 
metries are presented in Section [2j The notion of equivalence of nonclassical symmetries with 
respect to a transformation group or a set of admissible transformations plays a crucial role in 
our consideration and therefore is separately given in Section [3j Section O is devoted to review- 
ing known results on admissible transformations, point symmetries and equivalences in class (pQ), 
including discrete ones. The presentation of these results is important since they form a basis for 
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the application of our technique involving transformations between equations and are extended 
in the paper to both the classes of determining equations. Moreover, Lie symmetry operators 
are special cases of reduction operators. The determining equations are derived in Section [5] for 
both the cases of nonvanishing and vanishing coefficients of dt- It is proved in Section [6] via 
description of all possible reductions that solving the determining equations is equivalent to the 
construction of parametric families of solution of the corresponding initial equations. As a result, 
nonlocal transformations reducing the determining equations to the initial ones are found. Point 
transformations and Lie reductions of the determining equations are studied in Sections \7\ and El 
respectively. Results on Lie reductions of the determining equations corresponding to reduction 
operators with zero coefficients of dt are presented in such a form that they are directly extended 
to the general class of (1 + l)-dimensional evolution equations. In Section [9] we investigate the 
determining equations along with some non-Lie additional constraints. A nontrivial application 
of reduction operators to finding exact solutions of equations from class (pQ), arising under Lie 
reductions of the Navier-Stokes equations, is presented in Section [TUJ In the last section we dis- 
cuss possible extensions of obtained results, in particular, via study of the observed connection 
between reduction operators and the Darboux transformations of equations from class (PQ). 

To check results on Lie invariance of differential equations appearing in the paper, we used 
the unique program LIE by A. Head 

2 Reduction operators of differential equations 

Following [HI [35j [33] , in this section we shortly adduced necessary notions and results on non- 
classical (conditional) symmetries of differential equations. After substantiating with different 
arguments, we use the name "families of reduction operators" instead of "involutive families of 
nonclassical (conditional) symmetry operators" . 

Consider an involutive family Q = {Q 1 , . . . , Q 1 } of I (I ^ n) first order differential operators 

Q s = £ si (x, u)d l + 7] s (x, u)d u , 8 = 1,..., I, 

in the space of the variables x and u, satisfying the condition rank ||£ s *(x,u)|| = I. 

Hereafter x denote the n-tuple of independent variables (x±, . . . ,x n ) and u is treated as the 
unknown function. The index i runs from 1 to n, the indices s and a run from 1 to I, and we use 
the summation convention for repeated indices, di = d/dxi, d u = d/du. Subscripts of functions 
denote differentiation with respect to the corresponding variables. The local consideration is 
assumed. 

The requirement of involution for the family Q means that the commutator of any pair of 
operators from Q belongs to the span of Q over the ring of smooth functions of the variables x 
and u, i.e., 

Vs,s' 3C ssV = ( ss ' a (x,u): [Q s ,Q s '}=C ss ' a Q a . 

The set of such families will be denoted by Q l . 

If operators Q 1 , . . . , Q l form an involutive family Q, then the family Q of differential operators 

Q s = \ se7 Q a , where X S(T = X S(7 (x,u), det \\X sa \\ ^ 0, 

is also involutive and is called equivalent to the family Q. This will be denoted by Q = {Q s } ~ 
Q = {Q s }. (In the case I = 1 the functional matrix (A scr ) becomes a single nonvanishing 
multiplier A = \(x,u).) Denote also the result of factorization of Q l with respect to this 
equivalence relation by Ql. Elements of Q\ will be identified with their representatives in Q . 

If a family consists of a single operator (I = 1), the involution condition degenerates to an 
identity. Therefore, in this case we can omit the words "involutive family" and talk only about 
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operators. Thus, two differential operators are equivalent if they differ on a multiplier being a 
non- vanishing function of x and u. 

The first order differential function Q s [u] := r] s (x, u) — £ sl (x, u)U{ is called the characteristic of 
the operator Q s . In view of the Probenius theorem, the above involution condition is equivalent 
to that the characteristic system Q[u] =0 of PDEs Q s [u] = (called also the invariant surface 
condition) has n + 1 — I functionally independent integrals u>°(x,u), . . . ,uj n ~ l (x,u). Therefore, 
the general solution of this system can be implicitly presented in the form F(u)°, . . . , lo u ~ 1 ) = 0, 
where F is an arbitrary function of its arguments. 

The characteristic systems of equivalent families of operators have the same set of solutions. 
And vice versa, any family of n + 1 — I functionally independent functions of x and u is a 
complete set of integrals of the characteristic system of an involutive family of I differential 
operators. Therefore, there exists the one-to-one correspondence between Q\ and the set of 
families of n + 1 — I functionally independent functions of x and it, which is factorized with 
respect to the corresponding equivalence. (Two families of the same number of functionally 
independent functions of the same arguments are considered equivalent if any function from one 
of the families is functionally dependent on functions from the other family.) 

A function u = f(x) is called invariant with respect to the involutive operator family Q (or, 
briefly, Q-invariant) if it is a solution of the characteristic system Q[u] = 0. This notion is 
justified by the following facts. Any involutive family of I operators is equivalent to a basis 
Q = {Q s } of an /-dimensional (Abelian) Lie algebra g of vector fields in the space (x,u). Each 
solution it = f(x) of the associated characteristic system satisfies the characteristic system 
Q[u] = 0. Therefore, the graph of the function u = fix) is invariant with respect to the 
/-parametric local transformation group generated by the algebra q. 

Since rank ||£ sl (x, u)\\ = I, we can assume without loss of generality that ui^ / and F^o ^ 
and resolve the equation F = with respect to ll>°: lo° = tp{uJ l , . . . ,u n ~ l ). This representation 
of the function u is called an ansatz corresponding to the family Q. 

Consider an rth-order differential equation C of the form L{x,ur r \) = for the unknown 
function u of n independent variables x = (x\, . . . , x n ). Here, ur r ) denotes the set of all the 
derivatives of the function u with respect to x of order not greater than r, including u as the 
derivative of the zero order. Within the local approach the equation C is treated as an algebraic 
equation in the jet space of the order r and is identified with the manifold of its solutions 
in j( r ). Denote this manifold by the same symbol C and the manifold defined by the set of all 
the differential consequences of the characteristic system Q[u] = in by Q^ r \ i.e., 

Q (r) = {(x,u [r) ) e J (r) \D^ ...D^Q s [u] =0, oti G NU{0}, |a|: = «i + • • • + a n < r}, 

where Di = d Xi + u a+ s i d Ua is the operator of total differentiation with respect to the variable Xi, 
a = (ai, . . . ,a n ) is an arbitrary multi-index, 5i is the multiindex whose ith entry equals 1 and 
whose other entries are zero. The variable u a of the jet space corresponds to the derivative 
dWu/dx? 1 ...dx% n . 

Definition 1. The differential equation C is called conditionally invariant with respect to 
the involutive family Q if the relation Q s , r \L(x, u (r))\£ n Q(r) = holds, which is called the 
conditional invariance criterion. Then Q is called an involutive family of conditional sym- 
metry (or Q-conditional symmetry, nonclassical symmetry etc) operators of the equation C. 
Here the symbol Q 8 r r \ stands for the standard rth prolongation of the operator Q s [TTJ, l2Tj : 

Q s [r ) =Q S + E\aKr V sa du a , where V sa = • • ■ D^Q s [u] + $ si u a+Si . 

The equation C is conditionally invariant with respect to the family Q if and only if the 
ansatz constructed with this family reduces £ to a differential equation with n — I independent 
variables [44] . So, we will also call involutive families of conditional symmetry operators families 
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of reduction operators of C. Another treatment of conditional invariance is that the system 
C n is compatible in the sense of absence of nontrivial differential consequences [181 120J . 
Note that the paper [18] contains a number of interesting statements and ideas on the subject, 
which, unfortunately, did not become known according to their merits. If the infinitesimal 
invariance condition is not satisfied but nevertheless the equation C has Q-invariant solutions 
then Q is called a family of weak symmetry operators of the equation C [19, 20J. Nonclassical 
symmetries are often defined as generators of parametric groups of transformations preserving 
the solutions of C which additionally satisfy the corresponding invariant surface condition |12j . 
It is necessary to precisely interpret all the terms involved in this definition since otherwise it 
leads to the conclusion that, roughly speaking, any operator is a nonclassical symmetry of any 
partial differential equation. See also [HE] [20] f° r the discussion of connections between different 
kinds of symmetries. 

Lemma 1 ([Hllll]). If a differential equation is conditionally invariant with respect to an oper- 
ator family Q, then it is conditionally invariant with respect to any family of operators, which 
is equivalent to Q. 

The set of involutive families of I reduction operators of the equation £ is a subset of Q l and 
so will be denoted by Q l (C). In view of LemmaCD Q € Q l {£) and Q ~ Q imply Q € Q (£), i.e., 
Q l (C) is closed under the equivalence relation on Q l . Therefore, factorization of Q l with respect 
to this equivalence relation can be naturally restricted on Q l (C) that results in the subset Q\(C) 
of Ql. As in the whole set Ql, we identify elements of Q\{C) with their representatives in Q l (C). 
In this approach the problem of complete description of families of I reduction operators for the 
equation C is nothing but the problem of finding Q\{C). 

A different terminology can be used to call elements of Q\. Namely, it is possible to con- 
sider each element of Q\ as a C°°-module of the module dimension I, closed with respect to 
commutation [20], HI]. 

There are families of reduction operators related to classical Lie symmetries. Let q be an l- 
dimensional Lie invariance algebra of the equation C, whose basis operators satisfy the condition 
rank||£ SJ || = rank r/ s \\ ( = I' ^ I). The subsets consisting of I' elements of q, which are 
linearly independent over the ring of smooth functions of x and it, belong to Q l (£) and are 
equivalent each to other. The families of similar kind and ones equivalent to them will be called 
Lie families of reduction operators. The other families of reduction operators will be called 
non-Lie. 

3 Equivalence of families of reduction operators 
with respect to transformation groups 

We can essentially simplify and order the investigation of reduction operators, additionally 
taking into account Lie symmetry transformations in case of a single equation [25J and transfor- 
mations from the equivalence group or the whole set of admissible transformations in case of a 
class of equations [35] . Then the problem becomes similar to group classification of differential 
equations. 

Lemma 2. Any point transformation of x and u induces a one-to-one mapping of Q l into itself. 
Namely, the transformation g: x = X(x,u), u = U(x,u) generates the mapping gy. Q l — > Q l 
such that the involutive family Q is mapped to the involutive family g\Q consisting from the 
operators g*Q s = + fj s du, where £, st (x,u) = Q s X l (x,u), fj s (x,u) = Q s U(x,u). If Q' ~ Q 

then g\Q' ~ g\Q- Therefore, the corresponding factorized mapping g\: Q\ — > Q\ also is well- 
defined and one-to-one. 
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Definition 2 ( [32\ 125]). Involutive families Q and Q of the same number I of differential oper- 
ators are called equivalent with respect to a group G of point transformations if there exists a 
transformation g from G for which the families Q and g\Q are equivalent. 
Notation: Q ~ Q mod G. 

Lemma 3. Given any point transformation g of the equation C to an equation C, g\ maps Q l (C) 
to Q l (C) in a one-to-one manner. The same statement is true for the factorized mapping g\ 
from Q\{C) to Q\(C). 

Corollary 1. Let G be a Lie symmetry group of the equation C. Then the equivalence of involu- 
tive families of I differential operators with respect to the group G generates equivalence relations 
in Q l (C) and in Q\(C). 

Consider a class C\s of equations Cg: L(x,u^ r -j,9(x,u^)) = parameterized by 6. Here, L 
is a fixed function of x, ii( r ) and 9. The symbol 9 denotes the tuple of arbitrary (parametric) 
functions 9(x,u^) = (9 1 (x, «( r )), . . . , 9 k (x, U( r ))) running through the solution set S the sys- 
tem S(x, U( r ), U(r))) = 0. This system consists of differential equations on 9, where x and 
U( r ) play the role of independent variables and 9^ stands for the set of all the partial derivatives 
of 9 of order not greater than q. In what follows we call the functions 9 arbitrary elements. By 
G~ we denote the point transformations group preserving the form of the equations from C\s- 

For a fixed value / ^ n, consider the set P = P(L,S) of all pairs each of which consists 
of an equation Lq from L\s and a family Q from Q^{Cq). In view of Lemma El the action of 
transformations from G~ on C\s and {Q}{Cq) \ 9 £ S} together with the pure equivalence relation 
of involutive families of I differential operators naturally generates an equivalence relation on P. 

Definition 3. Let 9,9' £ S, Q £ Q l (£ s ), Q' £ Q l (C g >). The pairs (C S ,Q) and (C 9 ',Q') are 
called G~ - equivalent if there exists a transformation g £ G~ which maps the equation Co to the 
equation Cgi, and Q' ~ glQ- 

Classification of families of reduction operators with respect to G~ will be understood as clas- 
sification in P with respect to the above equivalence relation. This problem can be investigated 
in a way similar to the usual group classification in classes of differential equations. Namely, we 
construct first the reduction operators which are defined for all values of the arbitrary elements. 
Then we classify, with respect to the equivalence group, the values of arbitrary elements for 
which the corresponding equations admit additional families of reduction operators. 

In an analogous way we also can introduce equivalence relations on P, which are generated 
by either generalizations of usual equivalence groups or all admissible point transformations [30] 
(called also form-preserving ones [13]) in pairs of equations from £\$. 

Note 1. The consideration of the previous and this sections and known examples of studying 
reduction operators lead to the empiric conclusion that possessing a wide Lie symmetry group by 
a differential equation C complicates, in some way, finding nonclassical symmetries of C Indeed, 
any subalgebra of the corresponding maximal Lie invariance algebra, satisfying the transversality 
condition, generates a class of equivalent Lie families of reduction operators. A non-Lie family 
of reduction operators existing, the action of symmetry transformations on it results in a series 
of non-Lie families of reduction operators, which are inequivalent in the usual sense. Therefore, 
for any fixed value of I the system of determining equations on coefficients of operators from 
Q l (C) is not sufficiently overdetermined to be completely integrated in an easy way, even after 
factorized with respect to the equivalence relation in Q l (C). To produce essentially different 
non-Lie reductions, one have to exclude the solutions of determining equations, which give Lie 
families of reduction operators and non-Lie families being equivalent to others with respect to 
the Lie symmetry group of £. As a result, the ratio of efficiency of such reductions to expended 
efforts can be vanishingly small. 
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4 Lie group analysis of linear second-order parabolic equations 



Group classification in class (pQ) was first performed by S. Lie [2] as a part of his classification of 
general linear second-order PDEs in two independent variables. (See also a modern treatment 
of this subject in [21].) We shortly adduce these classical results, extending them for our pur- 
poses with using the notions of admissible transformations and normalized classes of differential 
equations. First, normalization properties of different classes of linear second-order parabolic 
equations were simultaneously analyzed in [34] in detail. 

Roughly speaking, an admissible transformation in a class of systems of differential equations 
is a point transformation connecting at least two systems from this class (in the sense that 
one system is transformed into the other by the transformation). The equivalence group of the 
class is the set of admissible transformations which can be applied to every system from the 
class. The class is called normalized if any admissible transformation in this class belongs to 
its equivalence group and is called strongly normalized if additionally the equivalence group is 
generated by transformations from the point symmetry groups of systems from the class. The 
set of admissible transformations of a semi-normalized class is generated by the transformations 
from the equivalence group of the whole class and the transformations from the point symmetry 
groups of initial or transformed systems. Strong semi-normalization is defined in the same 
way as strong normalization. Any normalized class is semi-normalized. Two systems from a 
semi-normalized class are transformed into one another by a point transformation iff they are 
equivalent with respect to the equivalence group of this class. See [251 [23 USB E2] for precise 
definitions and statements. 

Any point transformation T in the space of variables (t, x, u) has the form 

t = T t (t,x,u), x = T x (t,x,u), u = T u (t,x,u), 
where the Jacobian \d(T l , T x , T u )/d(t, x, u)\ does not vanish. 

Lemma 4. A point transformation T connects two equations from class (pQ) if and only if 
T x = T u = 0, T u = 0, T uu = 0, i.e., 

t = T(t), x = X(t,x), u = U 1 {t,x)u + U°(t,x), (2) 

where T, X, U 1 and U° are arbitrary smooth functions of their arguments such that TtX x U l ^ 
and additionally U° /U 1 is a solution of the initial equation. The arbitrary elements are trans- 
formed by the formulas 

T t T t \ U 1 J T t T t U 1 y J 

Here L = dt — Ad xx — Bd x — C is the second- order linear differential operator associated with 
the initial (non-tilde) equation. 

Corollary 2. Class JT]) is strongly semi-normalized. The equivalence group Gr~ of class ([1]) is 
formed by the transformations determined in the space of variables and arbitrary elements by 
formulas ([2]) and ([3]) ; where T, X and U 1 are arbitrary smooth functions of their arguments 
such that TtXxlI 1 ^ and U° = additionally. 

Note 2. Due to the presence of the linear superposition principle, the class ([I]) is not normalized 
because it is formed by linear homogeneous equations. The minimal normalized superclass of 
class ([I]) is the associated class of inhomogeneous equations of the general form 

u t = A(t, x)u xx + B(t, x)u x + C(t, x)u + D(t, x). 
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Using transformations from Gr~, the arbitrary elements A and B can be simultaneously 
gauged to 1 and 0, respectively. Hence, any equation from class (P) can be reduced by a 
transformation from G~ to an equation of the general form 

u t - u xx + V(t, x)u = 0. (4) 

The admissible transformations in the subclass are those admissible transformations in the 
class (pQ) which preserve the gauges A = 1 and B = 0, i.e., which additionally satisfy the 
conditions T/ = (T/ ) 2 and 22^Z* = -T t x T£. 

Corollary 3. A point transformation T connects two equations from class (jl]) if and only if it 
has the form 

t= fa 2 dt, x = ax + (, u = U 1 u + U°, U 1 := 9 exp f -^x 2 - — x ) , 

V 4<T 2a I 

(5) 

1 fir , aa tt ~ 2o- t 2 2 c< tt - 2a t ( t 8 t a t Q t 2 \ 

V A T* A T 



F = — —x 2 + 

a 2 V 4a 2 2a 2 2a 4a 2 J ' 

where a = o~(t), ( = C(t)> = an d = U°(t,x) are arbitrary smooth functions of their 
arguments such that aO ^ and U° /U 1 is a solution of the initial equation. Class is 
strongly semi-normalized. Any transformation from the equivalence group G~ of class (j3|) has 
the form ([5|) , where U° = additionally. 

The narrower equivalence group under preserving certain normalization properties suggests 
class flU as the most convenient one for group classification. Moreover, solving the group classi- 
fication problem for class ([1]) is reduced to solving the group classification problem for class (jlj) . 
The results on the group classification of class (pQ) (resp. (Jl|)) can be formulated in the form of 
the following theorem [14|, [2T] . 

Theorem 1. The kernel Lie algebra of class ([T]) (resp. is (ud u ). Any equation from class (pQ) 
(resp. @ ) is invariant with respect to the operators fd u , where the parameter-function f = 
f(t,x) runs through the solution set of this equation. All possible G~ -inequivalent (resp. G~ - 
inequivalent) cases of extension of the maximal Lie invariance algebra are exhausted by the 
following ones (the values ofV are given together with the corresponding maximal Lie invariance 
algebras): 

1. V = V{x): (d t , ud u , fd u ); 

2. V = fix~ 2 , fj,^0: (8 t , D, n, ud u , fdu); 

3. V = 0: (d t ,d x ,G, D,U,ud u , fd u ). 

Here D = 2td t + xd x , U = 4t 2 d t + Uxd x - (x 2 + 2t)ud u , G = 2td x - xud u - 

Let £ be an equation from class (pQ), g(C) denote its maximal Lie invariance algebra and 
q°°(C) be the infinite-dimensional ideal of this algebra, consisting of the operators of the form 
fd u , where the parameter-function / = f(t, x) runs through the solution set of C. The quotient 
algebra q(C)/q°°{C) is identified with the finite-dimensional subalgebra g css (C) of g(C), spanned 
by the 'essential' Lie invariance operators of C, which do not contain summands of the form 
f(t,x)d u - Each operator from g(C) is similar to an operator from g ess (C) under a trivial linear- 
superposition transformation t = t, x = x, u = u + f(t, x). 

Corollary 4. For every equation C from class ([T]) dim3 css (£) £ {1,2,4,6}. 

It will be shown below that for every equation £ from class ([I]) the number of reduction 
operators being inequivalent with respect to the group of linear-superposition transformations, 
roughly speaking, is significantly greater that the number of 'essential' Lie invariance operators. 
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5 Determining equations for reduction operators 
of linear second-order parabolic equations 

In the case of two independent variables t and x and one dependent variable u, each reduction 
operator is written as Q = r(t,x,u)d t + £(t, x, u)d x + n(t, x, u)d u , where (r, £) 7^ (0,0). The 
conditional invariance criterion for an equation £ from class (pQ) and the operator Q has the 
form [8] 

Q(2) Lu I l«=0, Q[u]=0, D t Q[u]=0, D x Q[u]=0 = °> 

where Q( 2 ) is the standard second prolongation of Q, Q[u] = n — riij — ^u x is the characteristic 
of Q, D t and D x denote the total differentiation operators with respect to t and x: 

D t = d t + u t d u + u u d ut + u tx d Ux H , 

As = d x + u x d u + u tx d Ut + i^aA^ H . 

All equalities hold true as algebraic relations in the second-order jet space 

J(2) 

over the space 

of the independent variables (t, x) and the dependent variable u. 

Since C is an evolution equation, there are two principally different cases of finding its re- 
duction operators: r ^ and r = 0. The investigation of these cases results in the preliminary 
description of the reduction operators. 

Lemma 5. Every reduction operator of an equation C from class ([TJ is equivalent to either an 
operator 

d t + g 1 (t,x)d x + (g 2 {t,x)u + g 3 (t,x))d u , 
where the functions g 1 = g 1 ^^), g 2 = g 2 (t,x) and g 3 = g 3 (t,x) satisfy the system 

g\ - Ag\ x - Bgl + Ugl - ^g 1 - ^jj) (g 1 + B) + B x g l + 2Ag 2 x + B t = 0, 

g\ - Ag 2 xx - Bgl + Ugl - ^-g 1 - ^) (g 2 - C) - C x g^ - C t = 0, (6) 

9! -A 9L -Bgl + (2 9 l-^-^y-C 9 3 = 0, 

or an operator d x + r/(i, x, u)d u , where the function r\ = r](t,x,u) is a solution of the equation 

r\ t = A(rj xx + 2rjrj xu + i] 2 i] uu ) + A x (rj x + rj7] u ) + (Bn) x + C(n - ur/ u ) + C x u. (7) 

Example 1. Each equation from class (pQ) with C = possesses the reduction operator d x . 

We denote the set of reduction operators of the equation £ from class JI]) by Q(C), omitting 
the superscript 1. The corresponding set factorized with respect to the equivalence of reduction 
operators is denoted by Qf (£). Consider the subsets Q\{C) and Qo(C) of Q(C), which consist of 
the operators constrained by the conditions r = 1 and (r, £) = (0, 1), respectively. The factor- 
set Q{{£) can be identified with Q\(C) U Qq(C). This union represents the canonical partition 
of <2f(£). The systems of form ([6]) and equations of form ([7]) associated with the equation C 
(and being the determining equations for the operators from Qi(C) and Qo(C)) are denoted by 
DEi(£) and DEo(£), respectively. It is obvious that the rules C — ► DEi(£) and C — > DEo(£) 
define one-to-one mappings of class (pQ) onto classes ([6]) and (J7]) . 
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Note 3. The partition of sets of reduction operators according to the condition of (non-)van- 
ishing of the coefficient r is natural for equations from class dU) (as well as the whole class of 
evolution equations) and agrees with their transformational properties. See Section [7] for details. 

Note 4. For certain reasons, here reduction operators are studied for equations of the non- 
reduced form (pTJ) . At the same time, it is enough, up to the equivalence relation generated by 
the equivalence group of class (pQ) on the set of pairs "(an equation of form ([1]), its reduction 
operator)", to investigate only the subclass (jl]) of equations with A = 1 and B = 0. The 
determining equations ([6]) and (J7]) for equations from class © have the simpler general form 

9 l t ~ 9 X XX + 2gW + ^al = 0, 

g 2 -g 2 xx + 2g 1 x (g 2 + V) + V x g 1 + V t = 0, (8) 
gl - gl x + 2glg* + Vg 3 = 

and 

Vt = Vxx + ZvVxu + V 2 Vuu ~ V(rj - ur] u ) - V x u. (9) 



6 Linearization of determining equations to initial ones 

There are connections between solution families of an equation C from class ([I]) and its reduction 
operators. This generates connections of the system DEi(£) and the equation DEo(£) with the 
initial equation C via nonlocal transformations. 

Consider at first reduction operators from Q\{C). Below the indices i and j run from 1 to 3. 
The indices p and q run from 1 to 2. The summation convention over repeated indices is used. 

Theorem 2. Up to the equivalences of operators and solution families, for any equation from the 
class (H|) there exists a one-to-one correspondence between its reduction operators with nonzero 
coefficients of dt and two-parametric families of its solutions of the form 

u = C\V 1 (t,x) +c 2 v 2 (t,x) + v 3 (t,x), (10) 

where c\ and C2 are constant parameters. Namely, each operator of such kind corresponds to 
the family of solutions which are invariant with respect to this operator. The problem of the 
construction of all two-parametric solution families of an equation ([]]), which are linear in pa- 
rameters, is completely equivalent to the problem of the exhaustive description of its reduction 
operators with nonzero coefficients of dt- 

Corollary 5. Non-linear coupled system (J6]) is reduced by the transformation 



19 19 1 9 I 9 

11 1) — 11 11 11 11 — 11 11 

1 _ a u u xx u xx u rj 2 A x xx u xx u x , /~t 

y ——/i — ; — 5 ; — 5 -D, g — — n. : — — 1- u, 
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A 



v l v 2 — v\v 2 



^ ^x ^xx 

v 2 v 2 v 2 

u w x u xx 

1 s s 

V v x v xx 



(11) 



to the uncoupled system of three copies of equation (fT]) for the functions v % = v l (t,x): 

Lv { = v\ - Av xx - Bv x - Ctf = 0, (12) 
and the functions v 1 and v 2 being linearly independent. 
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Note 5. Let W(ip l , . . . , ip n ) denote the Wronskian of the functions ip k = ip k (t, x), k = 1, . . . , re, 
with respect to the variable x, i.e. Wfy 1 , . . . , ip n ) = det(d l tp k /dx l )£ l=l . Then transforma- 
tion (jlip can be rewritten as 

!_ AW(v\v 2 )) x 2 _ Wjylvl) 3 _ W(v\v 2 ,v 3 ) 

9 W(v\v 2 ) 9 W{v\v 2 y ' 5 WV,« 2 ) ' 

The solutions = ip k (t,x), k = 1, . . . ,re, of an equation from class ([1]) are linearly independent 
if and only if W{tp l , . . . , (/j n ) ^ 0. See, e.g., Lemma 6 in [31j . Therefore, formulas (lllj) are well 
defined. 

Proof. Let £ be an equation from class ([!]) and Q = dt + g d x + (g 2 « + g 3 )<9« £ Qi(£), i.e., the 
coefficients g l = g l (t,x) satisfy the system DEi(£). An ansatz associated with Q has the form 
u = x)v9(<jj) + f°(t,x), where J 1 = f l {t,x) and /° = f°(t,x) are given coefficients, Z 1 7^ 0, 
(p = (p(io) is the new unknown function, cu = uj(t,x) is the invariant independent variable and 
<jj x 7^ 0. This ansatz reduces C to an (in general, inhomogeneous) linear second-order ordinary 
differential equation in 99, which we denote by C . The general solution of C is represented in 
the form 93 = c p ip p (io) + if 3 (to), where ip 3 is a particular solution of £ , 92 1 and ip 2 are linearly 
independent solutions of the corresponding homogeneous equation and c\ and C2 are arbitrary 
constants. Substituting the general solution of C! into the ansatz, we obtain the two-parametric 
family of solutions of C, having form (jlOp with v p = f(p p and v 5 = ftp 3 + g. The split in the 
equations Lu = and Q[u] =0 with respect to the constants c\ and C2 implies that each of the 
functions v l is a solution of C and 



{g l + B)vl-{g 2 -C)v p = -Avl 
{g x + B)v% - (g 2 - C)v 3 - g 3 = -An 



xx • 



Since v v 2 — v\v 2 7^ 0, the last system is a well-defined linear system of algebraic equations with 
respect to {g l ,g 2 ,g 3 ), whose solution is represented by (fTTl) . 

Conversely, suppose that J 7 is a two-parametric family of solutions of C, having form (llOj) . 
This means that each of the functions v l is a solution of C. The functions v 1 and v 2 are 
linearly independent since both the parameters c\ and C2 are essential. Consider the operator 
Q = dt + g d x + (g 2 u + g 3 )d u , where the coefficients g % are defined by (fTT|) . Q[u] =0 for any 
u E T . The ansatz u = v 1 ip{uj) + v 3 , where u> = v 2 /v , constructed with Q, reduces C to the 
equation <p uu = since (w 2 /^ 1 )^ = W (v 1 ,v 2 ) / {v 1 ) 2 / 0. Therefore 03j, Q G Qi(>C) and the 
functions 5* have to satisfy the system DEi(£). □ 

Corollary 6. Let C be an equation from class ([T]) and G°°(C) denote the trivial Lie invariance 
group of C, consisting of the linear superposition transformations of the form t = t, x = x and 
u = u + f(t,x), where the parameter-function f = f(t,x) runs through the solution set of C. 
Every reduction operator of the equation C with a nonvanishing coefficient of dt is G°°(C)- 
equivalent to an operator dt + g d x + g 2 ud u , where the functions g 1 = <7 1 (t, x) and g 2 = g 2 (t, x) 
satisfy the two first equations of DEi(C) . 

Proof. Suppose that a reduction operator Q of the equation £ has a nonvanishing coefficient 
of dt- In view of Lemma [5J the operator Q is equivalent to an operator Q of the form dt + 
g 1 d x + (g 2 u + g 3 )d u , where the functions g 1 = g l {t,x), g 2 = g 2 {t,x) and g 3 = g 3 (t,x) satisfy 
the system DEi(£). It follows from the proof of Theorem [2] that the coefficient g 3 possesses 
the representation g 3 = v 3 + g v 3 — g 2 v 3 , where v 3 = v 3 (t,x) is a solution of C. Then the 
transformation from G°°{C) with / = — v 3 maps the operator Q to the operator Q = dt + g 1 d x + 
(g 2 u + g 3 )du, where g 3 = g 3 — v 3 — g 1 v 3 + g 2 v 3 = 0. □ 
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Note 6. The functions v % satisfying the system f)12|) and the additional conditions (jlip with 
fixed values of the coefficients g 3 are defined up to the transformation 

v p = fi pq v q , v 3 = v 3 + fx 3q v q , (13) 

where fii q = const, and det(// P9 ) 7^ 0. Transformation (|13p induces the transformation of the 
constants c\ and C2: c p = fj- pq {c q — v q ), where (fi pq ) = (ii p > q i)~ l . It is obvious that the families of 
solutions (fTUj) and u = c~\v l + C2V 2 + v 3 coincides up to re-parametrization and can be identified. 

Consider reduction operators from Qq(C). 

Theorem 3. Up to the equivalences of operators and solution families, for any equation of 
form dH) there exists a one-to-one correspondence between one-parametric families of its solu- 
tions and reduction operators with zero coefficients of <9 t . Namely, each operator of such kind 
corresponds to the family of solutions which are invariant with respect to this operator. The prob- 
lems of the construction of all one-parametric solution families of equation ([1]) and the exhaustive 
description of its reduction operators with zero coefficients of c\ are completely equivalent. 

Corollary 7. The nonlinear (1 + 2) -dimensional equation (|7|) is reduced by composition of 
the nonlocal substitution rj = —& x /& u , where $ is a function of (t,x,u), and the hodograph 
transformation 

the new independent variables: t = t, x = x, x = 

the new dependent variable: u = u (14) 

to the initial equation Lu = in the function u = u(t, x, x) with x playing the role of a parameter. 

Proof. Let C be an equation from class ([]]) and Q = d x + r\d u G Qq(C), i.e., the coefficient 
77 = r](t,x,u) satisfies the equation DEo(£). An ansatz associated with Q has the form u = 
f(t,x,tp(uj)), where / = f l (t,x,Lp) is a given function, 7^ 0, 92 = ip(u>) is the new unknown 
function and lo = t is the invariant independent variable. This ansatz reduces £ to a first- 
order ordinary differential equation in <p, which we denote by £' . The general solution of CJ 
is represented in the form <p = (p(u,?t), where (p^ ^ and x is an arbitrary constant. The 
substitution of the general solution of £' into the ansatz results in the one-parametric family T 
of solutions u = f(t, x, x) of L with / = f(t, x, tp(t, x)). Expressing the parameter x from the 
equality u = f(t,x,x), we obtain that x = x,u), where $ M 7^ 0. Then t] = u x = —<& x /<& u 
for any u G i.e., for any admissible values of (t,x, x). This implies that r] = — <& x /& u for any 
admissible values of (t,x,u). 

Conversely, suppose that T = {u = f(t,x, x)} is a one-parametric family of solutions of C. 
The derivative f^ is nonzero since the parameter x is essential. We express x from the equality 
u = f(t,x,?c): x = &(t,x,u) for some function $ = $(t,x,u) with & u 7^ 0. Consider the 
operator Q = d x +rjd u , where the coefficient rj = r](t, x, u) is defined by the formula 77 = —$ x /$ u . 
Q[u] = for any u G T . The ansatz u = f(t, x, p(w)), where uj = t, associated with Q, reduces C 
to the equation <p w = 0. Therefore [44], Q G Qo(C) and hence the function 77 satisfies the equation 
DEo(£). □ 

Note 7. One-parametric families of solutions u = f(t, x, x) and u = f(t, x, x) of £ are assumed 
equivalent if they consist of the same solutions and differ only by parameterizations, i.e., if 
there exists a function ( = C,{x) such that ( x / and /(rj, x, £(x)) = f(t, x,x). Equivalent 
one-parametric families of solutions are associated with the same operator from Qq{£) and have 
to be identified. 

Note 8. The supposed triviality of the above ansatzes and reduced equations is connected with 
usage of the special representations for the solutions of the determining equations. Under this 
approach difficulties in the construction of ansatzes and the integration of reduced equations are 
replaced by difficulties in obtaining the representations for coefficients of reduction operators. 
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7 Admissible transformations, the equivalence groups 
and Lie symmetries of determining equations 



The "no-go" results of the previous section can be extended with investigation of point trans- 
formations, Lie symmetries and Lie reductions of determining equations (|6]) and (J7J). Thus, the 
maximal Lie invariance algebras of ([6]) and ([7]) are isomorphic to the maximal Lie invariance 
algebras of equation (pQ) in a canonical way. (Before this result was known only for the linear 
heat equation [7].) Moreover, the similar statements are true for the complete point symmetry 
groups including discrete symmetry transformations as well as the equivalence groups and sets 
of admissible transformations of classes of the above equations. 

All these statements are justified by Lemmas [3] and [H Indeed, each point transformation 
T between equations C and C from class ([!]) has the form ([2]) and induces the one-to-one 
mappings %: Q(C) — > Q{C) and 7J: Qf(£) — » Qf(£). Due to the conditions T x = and 
= 0, the transformation T* preserves the constraint r = (resp. r ^ 0) for coefficients 
of reduction operators. Therefore, the transformation 7f is split into the one-to-one mappings 
7f i: Qi(C) — ► Qi(£) and 7f : Qo(£) - ► Qo(£) according to the canonical partitions of Qf(C) 
and Qf(C). This implies that there exist the transformations T\ and To in the spaces of the 
variables (t, x, g 1 , g 2 , g 3 ) and (t, x, u, rj), which are induced by the transformation T in a canonical 
way. It is evident that 

Ti(DEi(£)) = DEi(£), %(DE (C)) = DE (£). 

The procedure of deriving the explicit formulas for 7{ is the following: Acting on the operator 
dt + g 1 d x + (g 2 u + g 3 )d u by T* and then normalizing the coefficient of &i to 1, we obtain the 
operator <9 t - + g l dx + (g 2 u + g s )du, where the new coefficients g l = g l (t,x), i = 1,2,3, are 
calculated by the formulas 

~i _ Xx_ i Xt 

9 ~ ~Y t 9 Tk' 



f = ±y + ^ + ^ d5) 



3 u 1 3 u° 2 u Q x u l -u Q u l x l u?u l -u*ut 

q = — Q or + — 1 — -q + — 1 — -■ 

Formulas (fTB]) describe the action of 7| on the dependent variables (g 1 , g 2 , g 3 ). The independent 
variables t and a; and the arbitrary elements A, B and C are transformed by the same formulas ([2]) 
and © as ones of the transformation T. The transformation of it is neglected. 

If the transformation T belongs to the equivalence group G~ of class ([I]) then it is defined 
for all values of arbitrary elements. Therefore, the same statement is true for 7i, i.e., T\ belongs 
to the equivalence group of class ([6]). In other words, the equivalence group of the initial 
class induces a subgroup of the equivalence group of the class of determining equations for the 
case t = 1. 

Suppose that the transformation T is parameterized by the parameter e and this family of 
transformations form a one-parametric Lie symmetry group of the equation C, generated by an 
operator Q = rdt + £,d x + (^u + (°)d u . We differentiate formulas (fT5l) with respect to e and 
then put e = 0, taking into account the conditions 

r = r(i) = TJ T\ r =t, £ = £(t,x) = X e \-, X\ = x, 



le=0' Ie=0 ' s svi"V ^£| e =0' le=0 

/^U, u\ = = i, c = c°(t,*) = ^u, u\ =Q - 
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As a result, we obtain the expressions for the coefficients 9 l of the Lie symmetry operator 
Qi = rdt + £,d x + 6 l d g i of the system DEi(£), associated with the operator Q: 

1 = (6* - r t )g l + i t , 

2 = -r t g 2 +rig l +r,l (16) 
e 3 = (r ? 1 -r,) 5 3 -r ? V+^ 1 +^ . 

The explicit formulas for Tq are derived in the analogous way. The action of on the operator 
d x + f]du and the normalization of the coefficient of d x to 1 result in the operator d x + rjd% , where 

_ U 1 U x U% 

T) = — V + ^u+-^. (17) 

A-x A-r A x 

Formula (|17p represents the expression for the dependent variable n transformed by 7o. The 
transformations of independent variables t, x and u and the arbitrary elements A, B and G are 
given by the formulas ((2J) and (jHJ). The unique difference from the transformation T is that the 
variable u is assumed independent. This implies that each transformation from the equivalence 
group G~ of class Q induces a transformation from the equivalence group Gq of class • 

Under the infinitesimal approach, each Lie invariance operator Q = rdt + £,d x + (C}u + C°)du 
of C is prolonged to the Lie invariance operator Qq = Q + 0d v of DEo(£), where the coefficient 9 
is determined by the formula 

e = (C 1 -Z*)v + & + & (is) 

The problem is to prove that the induced objects (resp. admissible transformations, point 
equivalences, point symmetries and Lie invariance operators) exhaust all possible objects of the 
corresponding kinds for determining equations. 

Lemma 6. If a point transformation connects two systems DEi(£) and DEi(£) from class (jSJ) 
then it has the form 

t = T{t), x = X(t,x), <f = G ii \t,x)g i ' + G i0 (t,x), (19) 

where T, X, G 33 and G 32 are smooth functions of their arguments such that TfX x G ss ^ 
and additionally G 32 /G 33 is a solution of the associated equation C; = 1,2,3. The other 
parameter-functions in (j!9p are explicitly defined: 

G 10 = ^, G" = ^, G 12 = 0, G 13 = 0, 

-^3F' G -jtfn> G ~T t > ° -°' (20) 

(T ^33 
f<30 _ \^t^ jt ^,31 _ ^ x ^,32 _ ^32 

T t 2 G 33 ' G 33 ' 

27ie arbitrary elements are transformed by the formulas 

A-Za, 6 = ^(b-^a)-^., C^-C^. (21 ) 

Here, L = dt — Ad xx — Bd x — C is the second-order linear differential operator associated with 
the equation C. 
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Proof. The systems DEi(£) and DEi(£) consists of second-order evolution equations which 
are linear in the derivatives, and coefficients of second derivatives form the nonsingular matrices 
diag(yl, A, A) and diag(A, A, A), respectively. In view of Corollary 13 of [M] each transformation 
between such systems necessarily has form f)19[) . We apply the direct method with taking into 
account the conditions (|19p and find more conditions which can be split by g % and g x . The 
system of determining equations on parameters of the transformation, obtained after the split, 
implies equations (|20l) and expressions ([3]) for transformations of the arbitrary elements. □ 

Theorem 4. There exists a canonical one-to-one correspondence between the sets of admissible 
transformations of classes ([ID and §5§. Namely, each point transformation between equations 
C and C from class ([T|) induces a point transformation between the associated systems DEi(£) 
and DEi(£) according to formulas (1151) . In both the transformations the independent variables 
are transformed in the same way. The induced transformations exhaust the sets of admissible 
transformation in class ©. 

Proof. It only remains to prove that every admissible transformation in class Q is induced by an 
admissible transformation in class (H|) in the above way. We fix two point-equivalent systems from 
class (J6j). They necessarily are systems of determining equations for reduction operators with the 
unit coefficients of dt for some equations C and C from class ([1]). Therefore, these systems can be 
denoted by DEi(£) and DEi(£), respectively. Consider a point transformation T mapping the 
system DEi(£) to the system DEi(£). In view of Lemma[6l the transformation T has form (|19p . 
where G 32 /G 33 is a solution of C and the other parameter- functions G u and G l ° are explicitly 
expressed by (|20p . Formulas (|2ip describe connections between the arbitrary elements of DEi(£) 
and DEi(£). We associate the transformation T with the transformation T in the space of the 
variables (t, x, u), having the form ([2]), where U 1 = T^G 33 and U° = TjG 32 . By the construction, 
U 1 /U is a solution of C. Since the pairs (DEi(£), DEi(£)) and (C, C) have the same tuples of ar- 
bitrary elements, LemmaHland formulas (|2ip imply that T is a point transformation from C to C. 
The comparison of (|20p with (|15p allows us to conclude that T is induced by T, i.e., T = T\. □ 

Note 9. It follows from the proof of Theorem U] that "if . . . then ..." in Lemma [6] can be 
replaced by ". . . if and only if . . .", i.e., the presented conditions are necessary and sufficient. 

Corollary 8. The equivalence group of class ([6]) is isomorphic to the equivalence group G~ 
of class (HJ). The canonical isomorphism is established by formulas (|15p . where U° = 0. 

Corollary 9. For each equation C from class ([T]), the maximal point symmetry groups (resp. 
the maximal Lie invariance algebras) of the equation C and the system DEi(£) are isomorphic. 
A Lie symmetry operator Q = rdt + £,d x + (C} u + C°)d u of C induces the the Lie symmetry 
operator Qi = rdt + £,d x + 9 l d g i of the system DEi(£), where the coefficients 9 l , i = 1,2, 3, are 
defined by formulas (|16p . 

Corollaries [8] and [9] along with Theorem [1] give the group classification of class © . 

Corollary 10. The kernel Lie algebra of class ([6]) is (1%), where I\ = g s d g 3. Any system from 
class Q is invariant with respect to the operators of the form Zi(f) = (ft + f x g l — fg 2 )d g 3, 
where the parameter-function f = f(t, x) runs through the solution set of the associated equation 
ft = Af xx + Bf x + Cf. All possible G^ -inequivalent cases of extension of the maximal Lie 
invariance algebra are exhausted by the following systems of the reduced form ([8]) (the values 
ofV are given together with the corresponding maximal Lie invariance algebras): 

1. V = V(x): (dt, h, Z 1 (f))- 

2. V = /ix" 2 , n + 0: (dt, Di, IT, h, Z 1 (f)); 

3. V = 0: (8 t , d x , G u Dt, Hi, I x , Z x (f)). 
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Here 



D 1 = 2td t + xd x - g'dgi - 2g 2 d g 2, 



Ui = 4t 2 d t + Uxd x + 4(x - tg l )d g i - (8tg 2 + 2xg 1 + 2)d g 2 - (x 2 + Wt)g 3 d g 
G x = 2td x + 2d Q i - g l d q 2 - xg 3 d Q 3. 



Note 10. It is obvious that Corollaries [U [9] and \W\ can be reformulated for subclass @ of 
the initial equations in the reduced form and subclass ([8]) of the corresponding determining 
equations of the first kind (the case r ^ 0) . 

A specific question for class © is what transformations of the functions (v , v 2 , v 3 ) defined in 
Corollary [5] are induced by admissible transformations in class ([6]) . It is clear that each induced 
transformation is admissible in class (|12p . Let £ and £ be equations from class (H|). Denote the 
corresponding systems of form (fT2l) by 3£ and 3£ and the corresponding sets of formulas (fTTI) 
by Q and Q, respectively. It is proved analogously to Lemma [6] that any point transformation 
connecting the systems 3£ and 3£ has the form 

t = T(t), x = X(t,x), v { = U 1 (t,x)fi ij v j + U i0 {t,x), 

where fMij = const, det(fiij) ^ 0, i, j = 1, 2, 3; T, X, U 1 and U t0 are arbitrary smooth functions of 
their arguments such that TtX x XJ x ^ and additionally U i0 /U 1 are solutions of the equation L. 
The arbitrary elements are transformed by the formulas (J3j) , where L = dt — Ad xx — Bd x — C 
is the second-order linear differential operator associated with the equation C. The agreement 
of transformations between 3C and 3C with transformations between DEi(£) and DEi(£) via 
formulas (jlip implies the additional conditions 

m=V23 = 0, U 10 = U 20 = 0, /x 33 = l, U^TtG 33 , U m =T t G 30 , 

for the admissible transformations between the systems 3£(~l^nDEi(£) and 3£n£?nDEi(£). To 
derive these conditions, we express all the tilde variables in Q via the non-tilde ones, then sub- 
stitute the expressions for g l given by Q into Q and split with respect to v 3 and their derivatives. 
Combining the obtained result with Theorem U] and omitting the systems DEi(£) and DEi(£) 
as differential consequences of the systems 3C(lQ and 3CnG, respectively, we get that the point 
transformation T of form ([2]) between the equations C and C induces the point transformation 

t = T{t), x = X(t,x), v p = U 1 (t,x)fi pq v q , v 3 = U 1 {t,x)fi 3q v q + U°(t,x), 

where det(^ pg ) ^ 0, p, q = 1,2, between the system 3C n Q and 3C n Q. The appearance of 
the additional constants fii q in the induced transformation is explained by the uncertainty (|13p 
under determining the function v % . The consideration of a one-parametric Lie symmetry group 
of the equation C instead of a single transformation between the (possibly different) equations 
C and £ results in a formula for the extension of Lie symmetry operators of C to Lie symmetry 
operators of 3C Namely, the following statement is true. 

Lemma 7. Each Lie symmetry operator Q = rdt+Cd x + {C lu ~ l rC )9u of the equation £ generates 
the family 

{rd t + id x + C l v l d v i + C,°d v 3 + \i q v q d v i \ \ iq = const} 

of Lie symmetry operators of the associated system 3£ with the additional conditions Q. Here 
i,j = 1,2,3, q = 1,2. The functions g 3 satisfy the system DEi(£) being the compatibility 
condition of 3£ H Q . 
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The chain of similar statements is also obtained for class (JT]). 

Lemma 8. // a point transformation in the space of the variables (t, x, u, rj) connects two equa- 
tions DEo(£) and DEo(£) from class ([7]) then it has the form given by formulas ([2]) and (fT?]) . 
where T, X, U 1 and U° are arbitrary smooth functions of their arguments such that TtX x U l ^ 
and additionally U° /U 1 is a solution of the equation C. The arbitrary elements are transformed 
by the formulas ([3]), where L = dt — Ad xx — Bd x — C is the second-order linear differential 
operator associated with the equation C 

Proof. The matrices formed by the coefficients of the second derivations in the equations DEo(£) 
and DEo(jC) are singular. That is why we cannot use the results of [36] on admissible transfor- 
mations in classes of parabolic equations having positively defined matrices of the coefficients of 
the second derivations. All determining equations have to be obtained independently. 

We use the direct method. Consider a point transformation T from the equation DEo(£) to 
the equation DEo(£) of the general form [t, x,u,fj] = [T, X, U, H] (t, x, u, rj) with the nonvanishing 
Jacobian. Sometimes we will also assume that the old variables (t, x, u, rj) are functions of the new 
variables (t, x, u, fj) and do a simultaneous split with respect to both the old and new variables. 
This trick is correct under certain conditions. We introduce the notations Q := D x + rjD u , 
Q := D x + fjD u and F := Qfj. In the old variables, the function F is expressed via t, x, u, rj, rj t , 
r\ x and rj u , and moreover {jF nt , F Vx , F nu ) / (0,0,0). (Indeed, the condition F nt = F Vx = F^ u = 
means that the function F depends only on (t, x, t, rj) in the old variables and, therefore, is a 
function of only (t, x, u, fj) in the new variables. Then we could split the equation F = fj x + fjfju 
defining F with respect to derivatives of fj and obtain the contradiction = 1.) 

The equation DEo(£) can be written in the form QF = . . . , where the right-hand side contains 
derivatives only up to order 1. We return to the old variables in DEo(£) and confine it to the 
manifold of the equation DEo(£), expressing the derivative rj xx from DEo(£) and substituting 
the found expression into DEo(£). Then we split the obtained equation DEq step by step with 
respect to different subsets of the other derivatives of n (or fj alternatively). To optimize this 
procedure, we start from the subsets of derivatives giving the simplest determining equations 
and take into account found equations for the further split. Note that the expression QF has 
the representation QF = (Qt)D t F + (Qx)D x F + (Qu)D u F. 

After collecting the coefficients of rju, rjt x and r/t u in DEq, we derive the system 

(Qt)F Vt = 0, (Qt)F Vx + (Qx)F Vt = 0, (Qt)F Vu + (Qu)F m = 

which implies the equation Qt = since (F Vt , F Vx , F^) ^ (0,0,0). We expand the expression 
Qt, assuming t a function of (t, x, u, fj): Qt = t x + tfjfj x + fjjtu + tfjfjxj). The split of the equation 
Qt = with respect to the new jet variables fj x and rju results in the equations tfj = and 
tx + fjta = 0. Then the subsequent split with respect to the new variable fj gives the equations 
t x = and t u = 0. Therefore, t is a function of only t, i.e., t depends only on t, t = T(t). Under 
this condition the function F expressed in the old variables does not depend on rjt, i.e., F^ = 
and hence (F rix ,F n J^ (0,0). 

Collecting the coefficients of ij uu and rj xu in DEq gives the system 

{Qu)F nu - v\Qx)F Vx = 0, (Qx)F Vu + (Qu)F Vx - 2n(Qx)F Vx = 0. 

Since (F^jF^) ^ (0, 0), the determinant of the matrix of this system considered as a system of 
linear algebraic equations with respect to (F Vx ,F Vu ) has to vanish, i.e., (Qu — tjQx) 2 = that 
implies Qu = r/Qx. Assuming x and u to be functions of (i, x, u, fj), we expand the expression Qx 
and Qu similarly to Qt and split the equation Qu = r\Qx with respect to the new jet variables 
fj x and rju- This results to the equations u^ = rjx^ and u x + fju u = rj(x x + fjxy). Alternating the 
old and new variables in any derived equation gives a correct equation. Therefore, we also have 
the equations U v = HX^, U x + rjU u = H(X X + rjX u ). 
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The next term for collecting coefficients in DEq is rj t . The equation obtained by this split is 
presented as AG = A(Qx)F rjx , where G denotes the coefficient of r/ t in fjj:. Under the above-stated 
conditions, the expressions appearing in this equation take the form 

1 ( D{H,U) u D(X,H) \ 1 8(H,X,U) ~ U u - HX U 

F = — — — — r + -" _ FT7 — r > = tt~a — o? r"> Q x ~ 



A\D(x,u) D(x,u) /' T t A d(r),x,u) ' A 

Hereafter A = D(X, U)/D(x,u) 0), and 
d(Z\...,Z fe ) ^ D(Z\...,Z k ) 



d(zi,...,Zk) D(zi,...,Zk) 

denote the usual and total Jacobians of the functions Z l , . . . , Z k with respect the variables z\, 
. . . , Zk, respectively. Note that in the case of a single dependent variable each total Jacobian 
is, at most, a first-order polynomial in derivatives of this dependent variable. Removing the 
denominators from the equations AG = A(Qx)F rix results in the equation 



o{r], x, u) 



A[ d{H 1 U) | H d(X,H) 
d(rj,u) d(r],u) 



d(X,U) fD(H,U) , H D(X,H) 



d(rj,u) \ D(x,u) D(x,u) 

the right-hand side of which is at most a first-order polynomial in rj x and r\ u . In view of 
nonvanishing A and d(H, X,U)/d(rj,x,u), this implies that the coefficients of r\ x and r\ u in A 
equal zero, i.e., d(X,U)/d(r],u) = and d(X,U)/d(x,r]) = 0. Then d(X,U)/d(x,u) / since 
otherwise the transformation T is singular. Hence X^ = 17^ = 0. 

Collecting the coefficients of r/ x in DE leads to the equation H m (U u — H X x ) 2 = 0. Note that 
U u — HX X = (Qx)A y£ since and Qx ^ 0. (Via the split with respect to unconstrained 

tilde variables, vanishing Qx implies the condition x x = x„ = Xfj = which contradict the 
nonsingularity of the inverse of T .) Therefore, H m = 0, i.e., H = H 1 ^, x, u)t] + H°(t, x,u), 
where H 1 = H v ^ 0. Knowing the explicit dependence of H on r\ allow us to additionally 
split all equations with respect to 77. Thus, splitting the equation U x + rjU u = H(X X + rjX u ) 
gives the condition X u = (hence X X U U 7^ 0) and, then, the conditions H 1 = U u /X x and 
H° = U x /X x . The equation DE contains only a single term including rj 2 rj u . Equating the 
corresponding coefficient to zero, we derive the condition U uu = 0. 

The whole set of the above found conditions on T, X, U and H implies that the form of the 
transformation 7" is described by formulas ([2]) and (|17p . Then the operator Q is transformed in 
a simple way: Q = X~ l Q. This gives us the idea to rewrite the equations DEo(£) and DEo(£) 
in terms of the operators Q and Q, respectively. Thus, the equation DEo(£) has the form 

rjt + r) u (AQrj + Br] + Cu) = AQ 2 r] + (A x + B)Qrj + (B x + C)r) + C x u. 

All derivatives of rj containing the differentiation with respect to x are excluded from DEq by 
the substitution rj x = Qu — r/r] u , and hence DEq can be split with respect to Q 2 r], r/ u , Qrj, rj 
and u. Collecting the coefficients of the terms r] u Qrj, r] u r], r\ u u and r] u , we obtain the formulas © 
for transformations of the arbitrary elements A, B and C and the condition L{U l /U°) = 0. □ 

Note 11. We do not split under deriving determining equations in the proof of Lemma [8] as 
much as possibly since the resulting system would be too cumbersome and, moreover, the proof 
of the next Theorem [5] implies that in fact this complete system is reduced to the set of conditions 
presented in Lemma El 
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Theorem 5. There exists a canonical one-to-one correspondence between the sets of admissible 
transformations of classes §\§ and ([7]). Namely, each point transformation between equations 
£ and £ from class ([T]) is extended to a point transformation between the associated equations 
DEo(£) andD~Eo(£) according to formula ()17p . In both the transformations the variables (t,x,u) 
and the arbitrary elements are transformed in the same way. The extended transformations 
exhaust the sets of admissible transformation in class (J7J). 

Proof. The extension of each admissible transformation in class ([1]) by formula (I17p gives an 
admissible transformation in class (JT]). Therefore, it is enough to check that every admissi- 
ble transformation in class ([7j) coincides with the extension of an admissible transformation in 
class (HJ). We take two equations from class © which are connected via a point transformation. 
They necessarily are determining equations for reduction operators with the zero coefficients of 
dt and the unit coefficients of d x for some equations £ and £ from class ([T]). Therefore, these 
equations can be denoted by DEo(£) and DEq(£), respectively. Consider a point transforma- 
tion T mapping DEo(£) to DEo(£). In view of Lemma [HJ the transformation T has the form 
given by formulas ([2]) and (fT7|) and, therefore, is projectable on the space of the variables (t, x, u). 
Denote its projection by T. The pairs (DEo(£), DEo(£)) and (£,£) have the same tuples of 
arbitrary elements transformed by the same formulas ([3]) . Hence Lemmas H] and [8] imply that 
T is a point transformation from £ to £. It is clear that the transformation T is the extension 
of T by formula (|17p. i.e., T = Tq. □ 

Corollary 11. The equivalence group Gq of class ([7]) is isomorphic to the equivalence group 
G~ of class ([T]) . The canonical isomorphism is established by the extension of transformations 
from Gq to the variable n via formula (|17p . where U° = 0. 

Corollary 12. For any equation £ from class ([I]), the maximal point symmetry groups (resp. the 
maximal Lie invariance algebras) of the equations £ and DEo(£) are isomorphic. The canonical 
isomorphism between the algebras is realized via the extension of each Lie symmetry operator 
Q = Tdt+£,d x + (( 1 u+(°)d u of £ to the Lie symmetry operator Q\ = Q+ ((C 1 — t.x)v+(x u +(x)dv 
o/DE (£). 

In view of Corollaries [TT] and [T2l the results on the group classification of class ([7]) follow 
from Theorem [TJ 

Corollary 13. The kernel Lie algebra of class ((7|) is {Io), where Iq = ud u + nd v . Any equation 
from class ([7]) is invariant with respect to the operators of the form Z${f) = fd u + fxd-q, where 
the parameter-function f = f(t,x) runs through the solution set of the associated equation 
ft = Af xx + Bf x + Cf . All possible Gq -inequivalent cases of extension of the maximal Lie 
invariance algebra are exhausted by the following equations of the reduced form Q (the values 
ofV are given together with the corresponding maximal Lie invariance algebras): 

1. V = V(x): I , Z {f)); 

2. V = /ix" 2 , v± 0: (dt, D , II , I , Z Q (f)); 

3. V = 0: (8 t , d x , Go, D , n , I , Z (/)>. 
Here 

D = 2td t + xd x - rjd v , 

n = At 2 d t + Atxd x - (x 2 + 2t)ud u - (xrj + 6tn + 2xu)d v , 
Gq = 2td x — xud u — (xn + u)d v . 
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8 Lie reductions of determining equations 



Suppose that an equation C from class ([I]) admits a Lie symmetry operator Q = rdt + £,d x + Qd u . 
The coefficients of Q necessarily satisfy the conditions t x = t u = 0, £ u = and ( uu = 0, i.e., 
r = r(i), £ = x) and C = C X (A x ) u + C°(*> x )' an d C° is a solution of £. 

In view of Corollaries l9l and 1121 the determining equations DEi(£) and DEo(£) respectively 
possess the Lie symmetry operators Q± and Qo associated with Q, which can be applied to 
reduce the determining equations and construct their exact solutions. The found solutions of 
the determining equations give the reduction operators of a special kind for the initial equation C, 
implicitly connected with Lie invariance properties of £ The question is what properties the 
solutions of C, invariant with respect to such reduction operators, possess, e.g., whether these 
solutions necessarily are Lie invariant or they are not. 

An admissible transformation T of the equation C in class ([I]) has form ([2]) and maps the pair 
(£, Q) to a pair (£,Q'), where the equation £ also belongs to class ([1]) and Q' is a nontrivial 
(resp. trivial) Lie symmetry operator of £ if Q is a nontrivial (resp. trivial) Lie symmetry 
operator of £ Up to the equivalence generated by the set of all admissible transformations of 
class dl]) (see Lemma H} in the set of pairs '(equation of form ([I]), its Lie symmetry operator)', 
we can assume that Q 6 {dt,d x } or Q E {ud u ,d u } if Q is a nontrivial or trivial Lie symmetry 
operator of £, respectively. Q ~ dt if r ^ and Q ~ d x if r = and £ ^ 0. 

If Q € {<9j, <9 X }, the Lie symmetry operator Qi of the system DEi(£) and the Lie symmetry 
operator Qq of the equation DEo(£), which are associated with the operator Q, formally have 
the same form as the operator Q but are defined in different spaces of variables. 

Proposition 1. Suppose that an equation C from class (pQ) possesses a Lie symmetry operator 
Q = rdt + £d x + (d u , where necessarily r = r(i), £ = £(i,x) and ( = C}{t,x)u + (°(t,x) and 
additionally r / 0. Let Q\ be the associated Lie symmetry operator of the system DEi(£), a 
solution G? 1 ,(? 2 ,g 3 ) o/DEi(£) be Qi-invariant and R = dt + g 1 d x + (g 2 u + g 3 )d u € Qi(£) be 
the corresponding reduction operator. Then the functions g 1 , g 2 and g 3 are expressed, according 
to formulas (fTTT) . via a solution (f 1 ,^ 2 ,?; 3 ) of the uncoupled system 3C, which is invariant with 
respect to the Lie symmetry operator 

rd t + £d x + C V^i + C 1 v 2 d v 2 + (C V + C°)^3 + X iq v q d vi 

of this system for some constants \i q , « = 1,2,3, q= 1, 2. Here the functions v 1 and v 2 have to 
be linearly independent. Each R-invariant solution of £ is a linear combination, with the unit 
coefficient ofv 3 , of the components of the Lie invariant solution (v 1 ,v 2 ,v 3 ) of the system 3£. 

Proof. It is sufficient to consider only the reduced form of Lie symmetry operators, which is 
Q = dt in the case r / 0. Then Qi = dt- The equation C is Lie invariant with respect to the 
operator dt if and only if At = B-t = Ct = 0. Consider an operator R = dt + g d x + (g 2 u + 
g s )d u G Qi(C), where the coefficient tuple (g 1 ,g 2 ,g 3 ) is a Qi-invariant solution of DEi(£), i.e., it 
additionally satisfies the condition g\ = g 2 = g\ = 0. An ansatz constructed with the operator R 
has the form u = f l (x)tp{uj) + f°(x), where f 1 = f 1 (x) ^ and f° = /°(x) are given coefficients, 
cp = tp(oj) is the new unknown function, u = t + g(x) is the invariant independent variable and 
g x ^ 0. This ansatz reduces C to an (in general, inhomogeneous) linear second-order constant- 
coefficient ordinary differential equation in <p, which we denote by £ . The general solution of £ 
is represented in the form ip = c p (p p (oj) + ip^ioo), where c^ 3 is a particular solution of £, ip 1 
and ip 2 are linearly independent solutions of the corresponding homogeneous equation and c\ 
and C2 are arbitrary constants. Let us recall that p,q = 1,2. Substituting the general solution 
of £ into the ansatz, we obtain the two-parametric family of solutions of C, having form (jlOD 
with v p = fip p and v 3 = fip 3 + g. Due to £ is a constant-coefficient equation, the functions v l 
admit the representation v p = tp pq (t)6 q (x) and w 3 = ip 3q (f)Q q {x) + 9 ?J (x), where ip iq = Xi p ip pq for 
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some constants Aj p depending on the coefficients of C Therefore (w 1 , v 2 , v 3 ) is a solution of the 
system 3£, which is invariant with respect to the Lie symmetry operator dt + \i q v q d v i of this 
system. □ 

Proposition 2. Suppose that the system DEi(£) associated with an equation C from class ([1]) 
possesses a Lie invariance operator Q\ with the vanishing coefficient of dt and a nonvanishing 
coefficient of d x . Let a solution (g 1 , g 2 , g 3 ) o/DEi(£) be invariant with respect to Q\. Then the 
associated reduction operator dt + g 1 d x + (g 2 u + g 3 )d u of the equation C is necessarily equivalent 
to a Lie invariance operator of C 

Proof. Consider the case Q = d x . The equation C possesses the Lie symmetry operator d x if 
and only if A x = B x = C x = 0. Then the equivalence transformation t = T(t), x = x + (p(t) 
and u = ip(t)u, where Tt = A, (ft = B ipt = Ctp and ip ^ 0, maps Q to d x and reduces C to the 
linear heat equation = u xx associated with the values A = 1 and B = C = 0. That is why 
without loss of generality we can assume that A = 1 and B = C = 0. An ansatz constructed 
for the system DEi(£) by the operator Q\ = d x is g % = g l {t) and the corresponding reduced 
system has the form g\ = 0, i.e., g % = const. The operator dt + g d x + (g 2 u + g s )d u with constant 
coefficients belongs to the maximal Lie invariance algebra of the equation C which coincides 
under our suppositions with the linear heat equation. The obtained statement is reformulated 
for the general form of Q with the vanishing coefficient of dt- □ 

Results on Lie solutions of the determining equation DEo(£) can be presented as a single 
statement without split into different cases depending on the structure of the corresponding Lie 
symmetry operators. To formulate them in a compact form, we need to introduce at first the 
auxiliary notion of one-parametric solution families of the equation £, associated with the Lie 
symmetry operator Q of C. The set of such families is partitioned into two subsets which are 
respectively formed by the singular associated families consisting of Q-invariant solutions of C 
and the regular associated families obtained via acting on fixed non-Q-invariant solutions of C 
by the one-parametric transformation group generated by Q. 

Let us recall that Qo denotes the Lie symmetry operator of DEo(£), associated with Q. 
Equivalent families of solutions, which differ only by parametrization, are identified. In particu- 
lar, regular one-parametric families associated with the same operator are equivalent if and only 
if they differ only by parameter shifts. Such families are obtained by the action of the same 
one-parametric transformation group on fixed solutions which are similar with respect to this 
group. A neighborhood of a nonsingular point of Q is considered. (Otherwise, the one-to-one 
correspondence in the next theorem may be broken. In some cases it can saved by taking into 
account discrete symmetry transformations, see Note 14 of |29].) 

Formulas (|17p and (|18p imply the following statement which will be used below. 

Proposition 3. Let an equation C from class ([T|) be invariant with respect to a point transfor- 
mation T (resp. an operator Q) and the function r] = rj(t,x,u) be a solution of the associated 
determining equation DEo(£). Then the equations u x = r)(t,x,u) admits the transformation T 
(resp. the operator Q) as a point symmetry transformation (resp. a Lie symmetry operator) if 
and only if the function rj is an invariant of the associated point symmetry transformation Tq 
(resp. the associated Lie symmetry operator Qo) of the equation DEo(£). 

Theorem 6. For each equation C from class ([1]) and each Lie symmetry operator Q of C, there 
exists a one-to-one correspondence between Qo-invariant solutions of the determining equation 
DEo(£) and one-parametric families of solutions of C, associated withQ. Namely, the reduction 
of the equation C by an operator d x + r)d u , where the coefficient r] is a Qo-invariant solution of 
DEo(£), gives a one-parametric solution family of C, associated with Q. And vice versa, each 
family of the above kind consists of solutions invariant with respect to an operator d x + rjd u , 
where the coefficient rj is a Q^-invariant solution of DEq(£) . 
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Proof. Suppose that an equation C from class (H|) admits a Lie symmetry operator Q. We 
denote the one-parametric transformation group with the infinitesimal operator Q by G. Let 
a solution n of the equation DEq(£) be invariant with respect to the associated operator Qq. 
Then the system of the equation C with the additional constraint u x = n possesses Q as a Lie 
symmetry operator. The general solution T of C v is a one-parametric solution family of C. There 
are two different cases of the structure of T. In the first case the family T consists of Q-invariant 
solutions of L and, therefore, is a singular one-parametric solution family associated with the 
operator Q. In the second case the family T contains a solution u = vP(t,x) of £, which is not 
Q-invariant. A one-parametric family of solutions of obtained via acting on the solution u° 
by transformations from G is equivalent to J-. Therefore, J- is a regular one-parametric solution 
family associated with the operator Q. 

Vice versa, if a one-parametric solution family of the equation C is associated with the 
operator Q then the corresponding additional constraint u x = n with a solution r/ of DEo(£) 
admits Q as a Lie symmetry operator. In view of Proposition [31 this implies that the function 
T) is Q-invariant. □ 

Since the determining equation DEo(£) has three independent variable, it also admits Lie 
reductions with respect to two-dimensional subalgebras of its maximal Lie invariance algebras to 
ordinary differential equations and, therefore, possesses the corresponding invariant solutions. 
To formulate the statement on such solutions analogously to Theorem El we need to define 
one-parametric families of solutions of the equation C, associated with the two-dimensional Lie 
invariance algebra g of C The whole set of associated families is also partitioned into the subsets 
of the singular and regular families. Each singular associated family consists of g-invariant 
solutions of C. Each regular associated family is obtained via acting on fixed (^-invariant and 
non-Q 2 -invariant solution of C by the one-parametric transformation group generated by Q 2 . 
Here Q 1 and Q 2 are arbitrary linearly independent elements of g. 

Theorem 7. Suppose that a two-dimensional Lie invariance algebra g of an equation C from 
class ([1]) induces the Lie invariance algebra Qq of the corresponding determining equation 
DEo(£), which is appropriate for Lie reduction o/DEo(£). Then there exists a one-to-one 
correspondence between go-invariant solutions of DEo(£) and one-parametric families of solu- 
tions of C, associated with g. Namely, the reduction of £ by an operator d x +rjd u , where the 
coefficient rj is a Qo-invariant solution of DEq(£), gives a one-parametric family of solutions 
of C, associated with q. And vice versa, each family of this kind consists of solutions invariant 
with respect to an operator d x + rjd u , where the coefficient n is a Qo-invariant solution o/DEo(£). 

Proof. We denote by G the two-parametric transformation group with the Lie algebra g and 
locally parameterize elements of G in a neighborhood of the identical transformation by the pair 
(ei,£2) : 5 f (£i>£2) £ G. In particular, g(0, 0) is the identical transformation and the infinitesimal 
operators Q 1 = g £i (0, 0), i = 1,2, form a basis of the algebra g. Let a solution rj of the equation 
DEo(£) be invariant with respect to the associated algebra go- Then g is a Lie invariance algebra 
of the system formed by the equation C and the additional constraint u x = rj. The general so- 
lution T of C v is a one-parametric solution family of C We explicitly represent this family by the 
formula u = f(t, x, x). There are two different cases of its possible structure. The family J- can 
consist of g-invariant solutions of C and, therefore, be a singular one-parametric solution family 
associated with the algebra g. The other possibility is that the family J- contains a solution 
u = f(t, x, hq) of which is not g-invariant. Then the solution u = f(t, x, *cq) is invariant with 
respect to the operator xo,\Q 2 — xo^Q 1 £ 0> where xo,j = (<?(ei, £2)x , o)eil(e 1 ,e 2 )=(o,o)- The action 
of the one-parametric subgroup G' of G with the infinitesimal operator xq^Q 1 + xq^Q 2 G g is 
(locally) transitive on T. It means that T is a regular one-parametric solution family associated 
with the algebra g, which is obtained via acting by G' on the fixed solution u = f(t,x,xo). 
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Conversely, if a one-parametric solution family of the equation C is associated with the 
algebra then the corresponding additional constraint u x = rj, where r] = rj(t, x, u) is a solution 
of DEo(£), admits g as a Lie symmetry algebra. In view of Proposition [3l this implies that the 
function rj is g-invariant. □ 



9 Particular cases of reductions and linearization 

All possible Lie reductions of the determining equations ([6]) and ([7]) for reduction operators and 
the corresponding invariant solutions of equations from class ([1]) are studied in the previous 
section. Now we consider a few examples of typical additional conditions to the determining 
equations, which are different from Lie ones. A special attention is paid to an interpretation of 
the confinement of the linearizing transformations given in Corollaries [5] and [7] to the particular 
cases under consideration. Presented examples also show that nontrivial reduction operators 
associated with nontrivial additional conditions to the determining equations can finally leads 
to trivial solutions of equations from class ([1]). 

We fix an equation C from class ([T|). The extension of possibilities for constraints of the 
determining equations in comparison with the initial equation C is connected with a greater 
number of unknown functions in DEi(£) and the additional independent variable u in DEo(£). 

Consider at first reduction operators of C with the vanishing coefficients of dt- 

Example 2. Suppose that Qq = d x is a reduction operator of DEo(£). It means that the arbi- 
trary elements satisfy the condition A x = B xx = C x = 0. The problem is to investigate solutions 
of DEo(jC), which are invariant with respect to Qq. We do an equivalence transformation of the 
form 

t = T{t), x = X 1 {t)x + X°(t), u = U l {t)u, 

where the arbitrary elements A, B and C and the function rj are transformed according to 
formulas ([3]) and (|17p . The parameter- functions T, X , X° and U 1 can be chosen in such a 
way that A = 1, B = and C = 0. In the new variables the operator Qq equals X d x and 
hence is equivalent to d x . This is why we can assume without loss of generality that A = 1, 
B = and C = 0, i.e., C coincides with the linear heat equation. Then Qq = d x is a Lie 
symmetry operator of DEo(£). The corresponding reduced equation rjt = r\r\ uu for the function 
rj = rj(t, u) is equivalent, on the subset of nonvanishing solutions, to the remarkable nonlinear 
diffusion equation Q t = (C~ 2 (u)u, where ( = 1/rj. It is well known that this diffusion equation 
is linearized to the linear heat equation [H [39] . We derive this transformation via confining the 
transformation of DEq(£) to, formally, C, presented in Corollary [71 We put $ = u) — x, 
where ^ 0. Then 

- J_ 

i.e., C = ^u- After integrating, we obtain the equation fy t = ^ uu /^ u 2 + f3(t) in the function 
$ = *f?(t, u). The "integration constant" j3 = (3(t) can be assumed to vanish due to the ambiguity 
in the connection between £ and The confinement of transformation (|14p is the hodograph 
transformation 

the new independent variables: t = t, x = ^, 
the new dependent variable: u = u 

since here the variable x has to be replaced by ^ = x + The application of this transfor- 
mation results in the linear heat equation = u xx . Note that the above interpretation of the 
confinement of transformation (|14p differs from the interpretation in the proof of Theorem 9 
of [29]. 
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Example 3. Let the function n satisfy the additional condition r\ uu = 0, i.e., r\ = r] 1 (t,x)u + 
r]°(t,x). Then the equation DEo(£) is reduced to the system 

T^ = {Ar^+A(r, 1 ) 2 + Bri L + C) w , 

rg = A( V ° XX + 2-rP-nl) + Mnl + A 1 ) + (Br] ), + Cri°. 

Putting = x)u-\-&°(t, x), we rewrite the transformation described in Corollary [7] in terms 
of r] 1 and rf . The condition rj = —& x /& u implies that rj 1 = — t&i/tfr 1 and rp = — The 
hodograph transformation (|14p is equivalent to expressing u from the formula for <3?: 

where V? 1 = l/^ 1 and ^° = $°/$ . Since the expression for it has to be the solution family 
of £ with the parameter x = $, ^ and \E , ° are solutions of C, ^f 1 7^ 0. Finally we derive the 
representation 

^ = |f, 17° = ^-^*°, (23) 

where $ 1 and \I /0 are solutions of the initial equation C. In other words, transformation (|23|) 
reduces the nonlinear system (|22[) in r/ 1 and 7/ to the system of two uncoupled copies of C. 
The expression for n 1 in f|23[) coincides, up to sign, with the well-known Cole-Hopf substitution 
linearizing the Burgers equation. (If A = 1 and B = C = 0, the first equation of (|22H coincides, 
up to signs, with the Burgers equation.) The expression for n 1 in (|23p is obtained as the Darboux 
transformation of the solution \& , associated with the solution . It follows from (123p that 
the reduction operator R = d x + (77 u + r/°)<9 u is G°°(>C)-equivalent to the operator d x + r] l ud u . 
Indeed, the transformation t = t,x = x,u = u — *° belongs to G°°(£) and maps the operator R 
to the operator R = d x + rj ud%. 

An ansatz constructed with i? has the form u = x)<p(uj) + x), where 99 = 95(0;) 
is an invariant unknown function of the invariant independent variable u = t. The associated 
reduced equation is ip^ = 0, i.e., tp = const. Therefore, u = ^ ,1 x+ l I /0 is the family of i?-invariant 
solutions of C. 

Vice versa, the solution family u = \I' 1 (t, x)n+ ^°(t, x) of the equation C is necessarily invari- 
ant with respect to the reduction operator d x + (r] 1 (t, x)u + r]° (t, x))d u , where the coefficients rj 1 
and rp are determined by the formulas (|23f) . 

As a result, we obtain the following statement. 

Proposition 4. For any equation of form (pQ), there exists a one-to-one correspondence between 
one-parametric families of its solutions, linearly depending on parameters, and reduction opera- 
tors of the form d x + r](t,x,u)d u , where r] uu = 0. Namely, each operator of such kind corresponds 
to the family of solutions which are invariant with respect to this operator. 

Example 4. At first sight, the additional condition rj x + rjr] u = seems much more complicated 
than the conditions studied in the previous examples. In fact, it leads only to solutions of the 
initial equation C, which are first-order polynomials with respect to x. To see this, we carry out 
the transformation described in Corollary [7J and, as a result, obtain the condition u xx = 0. In 
contrast to the solutions of C, the associated solutions of DEo(£) have a complex structure and 
are difficult to construct. 

The system S consisting of the equations DEo(£) and r\ x + r]7] u = has the compatibility 
condition (B xx + 2C x )r\ + C xx u = 0. Before considering the possible cases, we note that the 
equation rj x +rjrj u = is invariant with respect to the transformations from the equivalence group 
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Gq of class ©, which additionally satisfy the conditions {Ul/{U l ) 2 ) x = and {X x /{U l ) 2 ) x = 0. 
Denote the subgroup of Gq, consisting of these transformations, by Gq . The solutions of the 
system S are constructed up to G^-equivalence. 

If B xx + 2C X ^ 0, the function r\ has the form r\ = 77 1 (t, x)u. Then rj 1 = and C x = up to 
Gq -equivalence. The interpretation of this solution is obvious. An associated ansatz for C and 
the corresponding reduced equation are u = (f(u)), where to = t, and = 0. The family of the 
associated invariant solutions of C is formed by the constant functions. 

The condition B xx + 2C X = implies C xx = 0. Up to Gq -equivalence we can assume that 
B = C = 0. Then the system S is reduced to the system r\t = 0, rj x + r]rj u = 0. Its nonzero 
solutions are implicitly determined by the formula u = xrj + w(rj), where w = 10(77) * s an 
arbitrary smooth function of 77. An associated ansatz for the equation C is found from the 
condition u = xu x + w(u x ) which is the Clairaut's equation with the implicit parameter t. We 
choose the ansatz u = (p(u>)x + w(ip(u>)), where to = t. The corresponding reduced equation is 
<*Pu = 0, i.e., the associated invariant solutions of C has the form u = cx + w(c), where c is an 
arbitrary constant. 

Let us emphasize that the obtained results have a compact form only due to the consideration 
up to Gq -equivalence. 

Now we present a single example concerning the system DEi(£). In view of Corollary [6] we 
can assume without loss of generality that g 3 = and, therefore, consider only the two first 
equations of the system DEi(£). The G^-invariance of the equation g 3 = additionally justifies 
this assumption. 

Example 5. The constraint g 2 = is invariant with respect to the transformations from the 
equivalence group G~, in which U 1 = 1. These transformations are presented by formulas ([2]), 
([3]) and (fT5|) . where U 1 = 1 and U° = 0, and form the subgroup of G±, denoted by G±. Up 
to the G5"-equivalence, the coefficient A can be assumed equal to 1. Imposing the conditions 
g 2 = g 3 = and A = 1, we reduce DEi(£) to the system 

g\ - g\ x + 2g l g l x + (Bg l ) x + B t = 0, (24) 

C t +g l C x + 2g 1 x C = 0. (25) 

Equation (|24[) is linearized to the equation wt = w xx + (Bw) x by the generalization g 1 = 
—w x /w — B of the Cole-Hopf substitution and then to the equation vt = v xx + Bv x by the 
subsequent substitution w = v x . In the case C = 0, the resulting substitution g 1 = —v xx /v x — B 
is the confinement of transformation (jllj) under the assumptions v 3 = 0, v 2 = 1 and v 1 = v , 
where f is a nonconstant solution of C. 

Equation ([25]) admits a double interpretation depending on a reading of the phrase "the 
equation C possesses the reduction operator dt + g d x " . It can be considered either as an 
additional constraint for the function g 1 or an equation in the coefficient C. Choosing the 
second alternative, we obtain C = v x <&{v) for some function $ = 3>(u). 

If C = 0, equation ([25]) is an identity. Therefore, the equation C admits any reduction 
operators of the form dt — (v xx /v x + B)d x , where v = v(t, x) runs through the set of nonconstant 
solutions of C. The corresponding two-parametric solution family of C is u = c±v(t,x) + C2. 

Note 12. Since we do not initially specify values of the arbitrary elements and derive conditions 
on arbitrary elements depending on possessed reduction operators, the above examples have fea- 
tures of inverse problems of group analysis. Namely, we simultaneously describe both reduction 
operators with certain properties and values of arbitrary elements for which the corresponding 
equations admit such reduction operators. A similar inverse problem for generalized conditional 
symmetries of evolution equations is investigated in [38J. Due to possibilities on variation of 
arbitrary elements and application of equivalence transformations, the problems of this kind 
essentially differ from the problem of finding reduction operators of a fixed equation. 
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10 Applications 



In Sections [6HS] "no-go" statements of different kinds have been proved for the reduction opera- 
tors of the equations from class dTJ. The term "no-go" has to be treated only as the impossibility 
of exhaustive solving of the problem or the inefficiency of finding Lie symmetries and Lie reduc- 
tions of the determining equations. At the same time, imposing additional (non-Lie) constraints 
on coefficients of reduction operators, one can construct particular examples of reduction oper- 
ators and then apply them to the construction of exact solutions of an initial equation. Since 
the determining equations have more dependent or independent variables and, therefore, more 
degrees of freedom than the initial ones, it is more convenient often to guess a simple solution or 
a simple ansatz for the determining equations, which can give a parametric set of complicated 
solutions of the initial equations. (A similar situation is for Lie symmetries of first-order ordi- 
nary differential equations.) It is the approach that was used, e.g., in [TO] to construct exact 
solutions of a (nonlinear) fast diffusion equation with reduction operators having the zero coeffi- 
cients of dt- Earlier this approach was applied to the interesting subclass of class ([I]), consisting 
of the linear transfer equations of the general form 

h(t) , , 

ut = u xx H u x . (26) 

x 

These equations arise, in particular, under symmetry reduction of the Navier-Stokes equations 
[5j [22l E3] • Investigation of reduction operators allowed us to construct series of multi-parametric 
solutions of equations (|26p and, as a result, wide solution families of the Navier-Stokes equations, 
parameterized by constants and functions of t. 

We consider class (I26p as an example showing possible ways of imposing nontrivial additional 
constraints to determining equations. This subclass is singled out from the whole class ([1]) by 
the conditions on arbitrary elements A = 1, (xB) x = and C = 0. 

We fix an equation C from class (|26|) . The maximal Lie invariance algebra of C is the algebra 

1) (ud u , fd u ) if h / const; 

2) (d t , D, U h , ud u , fd u ) ifh = const, h {0,2}; 

3) {dt, D, U h , 28 x - hx~ l ud u , G h , udu, fdu) if h € {0,2}. 

Here D = 2td t + xd x , U h = At 2 d t + 4txd x - (x 2 + 2(1 + h)t)ud u , G h = 2td x - (x + htx'^udu. 
The function / = f(t, x) runs through the set of solutions of C. The case h = 2 is reduced to the 
linear heat equation (h = 0) by the transformation t = t, x = x and u = xu, cf. Theorem[TJ The 
inversection of the maximal Lie invariance algebras of equations from class (126)) coincides with 
(ud u , d u ), i.e., the kernel Lie symmetry group of class (|26p consists of scalings and translations 
of u. 

It is easy to see that the equation C possesses no nontrivial Lie symmetries and, therefore, 
no Lie reductions if h ^ const. At the same time, non-Lie reduction operators can be found for 
an arbitrary value of h. 

Any reduction operator of C with the nonzero coefficient of dt is G°°(£)-equivalent to an 
operator dt + g 1 d x + g 2 ud u , where the functions g l = g (t,x) and g 2 = g 2 (t,x) satisfy the 
two first equations of the corresponding determining system DEi(£). Following Example El we 
impose the additional constraint g 2 = 0. Then the second equation of DEi(£) is identically 
satisfied. The first equation of DEi (£) is rewritten in the form 

(g 1 + hx- l ) t = [gl - g 1 {g 1 + h x - 1 )) x . 

We put the left and right hand sides equal to 0. Theng 1 = x( x )~ hx~ l , g\— g (g -\-hx _1 ) = ip(t). 
The compatibility of these equations implies that x = —x -1 and ■0 = 0, i.e., g 1 = — (h(t) + l)x _1 
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and the corresponding reduction operator is 

Q = d t - (h(t) + 1^ 1 ^. 

As a result, the equation C possesses the family of Q-invariant solutions 

u = c 2 (x 2 + 2 J(h(t) + l)dt \ + ci. (27) 

Each reduction operator of C with the zero coefficient of dt is equivalent to an operator 
d x + f]d u , where the coefficient r] = rj(t,x,u) satisfies the corresponding determining equation 
DEoOC): 

Vt = Vxx + ZriVxu + r] 2 Vuu + h(x~ l rj) x . (28) 

Suppose that the same operator d x + r]d u is a reduction operator of all equations from class (i26j) . 
i.e., the function is a solution of (|28p for any value of h. This demand leads to the additional 
constraint (x~ 1 t]) x = implying that rj = x((t,u). We substitute the expression for rj into (|28l) 
and split with respect to x. Integrating the obtained system ( uu = 0, (t = 2CCu> we construct 
all its solutions: 

£ = ~^77^ — \ or C = 
2{t + x) 

where [/,, x and are arbitrary constants. In other words, the common reduction operators 
of equations from class ([26]) are exhausted, up to equivalence with respect to the kernel Lie 
symmetry group (more precisely, up to translations of u), by the operators of the form 

Gx = (2t + x)d x — xud u and d x + vd u . 

(It is obvious that there are no common reduction operators with nonzero coefficients of dt-) 
The constant x cannot be put equal to similarly to the constant fi since translations of t do not 
belong to the kernel Lie symmetry group of class (f26"J) and the classification up to the equivalence 
group of class (|26p in not convenient for the consideration. The operator is represented as 
the linear combination G + >cd x of the Galilean operator G = 2td x — xud u and translational 
operator d x . The non-reduced form for the coefficient of d x in G^ is chosen to obtain this 
representation. For any equation C from class (|2"6j) the reduction operator R = d x + ud u is 
G°°(jC)-equivalent to the operator d x which is trivial since the arbitrary element C equals in 
class (126(1 . Another formulation of above result is the following: Each equation from class (|26|) 
is conditionally invariant with respect to arbitrary linear combinations of the Galilean operator 
G and the translational operator d x . The family of G^- invariant solutions of an equation of the 
form (|26p consists of the functions 

f x 2 [ Ht) + l u 

u = ci exp < ; r — / at 

P \ 2(2i + x) J 2t + x 

The corresponding family for the operator d x + vd u has the form (|27p with c 2 = v. 
The constructed exact solutions are generalized to series of similar solutions 

k=0 fc=0 \ 1 / ' 

The functions T k = T k (t) and S k = S k (t) respectively satisfy systems of ODEs 

T t k = 2(k + l)(h(t) + 2k + l)T k+1 , k = 0,N -1, Jf = 0, 

= 2(k + l){h(t) +2k + l){2t + x)~ 2 S k+1 , k = 0,N-l, 5f = 0, 

which are easily integrated. These series of exact solutions also can be found using different tech- 
niques connected with reduction operators and their generalizations, in particular, via nonlocal 
transformations in class (I26p . associated with reduction operators [5| I22j. 
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11 Discussion 



The main result of the present paper is the chain of "no-go" statements on reduction opera- 
tors of linear (1 + l)-dimensional parabolic equations. These statements shows that application 
of conventional methods to solving of the determining equations for coefficients of such oper- 
ators cannot lead to reduction operators giving new exact solutions of initial equations. In 
both the cases naturally arising under the consideration, the determining equations form well- 
determined systems whose solving is in fact equivalent to solving of the corresponding equations 
from class ((TJ. All transformational and symmetry properties of the determining equations 
are induced by the corresponding properties of the initial equations. Reduction operators con- 
structed via Lie reductions of the determining equations are also connected with Lie invariance 
properties of the initial equations. Nevertheless, it is demonstrated in Section [10] that involve- 
ment of ingenious empiric approaches different from the Lie one can give reduction operators 
which are useful for the construction of non-Lie exact solutions of equations from class ([I]). 

Techniques developed in this paper can be applied to the general class of (1 + l)-dimensional 
evolution equations. We also plan to consider generalized reduction operators of linear (1 + 1)- 
dimensional parabolic equations, whose coefficients depend on derivatives of u. An interesting 
subject related to this is the connection between (generalized) reduction operators and Darboux 
transformations. Here we give some hints on this connection. 

Consider a fixed tuple of linearly independent functions (ip 1 , . . . , ip p ) of t and x, and the linear 
independence are assumed over the ring of smooth functions of t. The Darboux transformation 
constructed with the tuple (ip , . . . , ip p ) is denoted by DT[ifi , . . . ,ip p ] and is defined by the 
formula UHl 



u = DTW,...,^](u)= WW ^ ¥) • 

Here W((p l , . . . , (p s ) denote the Wronskian of the functions </? , <p s with respect to the 
variable x, i.e., W{tp l , . . . , Lp s ) = det{d % ~ 1 ^ j '<9x l_1 )jj =1 . The initial (u) and, therefore, obtained 
(u) functions also depend on t and x. 

The transformation DT^ 1 , • • • , ip p ] is represented as the action of a linear p-order differential 
operator with differentiations with respect to only x, DT^ 1 , . . . , ip p \(u) = DTf^ 1 , . . . , ip p \u. 
The operator will be denoted by the same symbol as the transformation and called the Darboux 
operator associated with the tuple (ip , ■ ■ ■ , ip p ). In the cases p = 1 and p = 2 the expressions of 
the Darboux operators respectively are 



If the functions ip 1 , . . . , ip p are linearly independent solutions of an equation C from class ([1]) 
then they are linearly independent over the ring of smooth functions of t |31|, I34j. The Darboux 
transformation DTf^ 1 , . . . ,ip p ] maps the equation C to the equation C also belonging to the 
class dl]) and having the following values of arbitrary elements |16j [M] 



A = A, B = B + P A X , C = C + pB x + p{p + 1} A xx + ^A x + 2 {^j A, 

where the abbreviation W = W(ip , . . . , ip p ) is used. 

Suppose that a reduction operator Q of L has the canonical form and is associated with a 
first-order linear differential operator Q acting on functions of t and x. It means that either 
Q = d t + g 1 d x + g 2 ud u if Q G Qi(£) or Q = d x + rfudu if Q G Qo(£)- (Here g 1 , g 2 and n 1 are 
functions of t and x.) In the first case the operator Q = —dt — g d x + g 2 equals the operator 
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— ^4 DTfu 1 , t> 2 ] on the solution set of the equation £, where the solutions v i = v i (t,x), % = 1,2, 
of C are determined according to Corollary [5j In the second case the coefficient n 1 admits the 
representation r/ 1 = ^f x /^>, where \l/ = ty(t,x) is a solution of C. Therefore, Q = — DT[^]. 
Finally, we have the following statement. 

Proposition 5. Let a reduction operator Q of an equation C from class (pQ) be associated, up 
to the equivalence relations of operators, with a first-order linear differential operator acting on 
functions oft and x. Then it is equivalent to a Darboux operator constructed with one (resp. two) 
linearly independent solutions of this equations in the case of vanishing (resp. nonvanishing) 
coefficient of dt- 

The properties of single reduction operators of multi-dimensional equations essentially differ 
from that in the (1 + l)-dimensional case. For example, all single reduction operators of (1 + n)- 
dimensional linear heat equations are exhaustively classified in [32J for arbitrary n without 
addressing the general solution of this equation that annuls the possibility of "no-go" statements. 
At the same time, it is not the case for involutive families of reduction operators [24| 140] . 
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